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Abstract 

The pressure - temperature phase diagrams of Se and S are calculated. 
Both melting and polymorphous phase transitions are described in frames 
of statistical Patashinski model. The results are in good agreement with 
experimental data of Brazhkin et. al. 
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1 Inroduction. 

Structural transitions in liquid Se and S similar to those in crystalline state 
were discovered by Brazhkin et.al. |^]. At these transitions, drastic changes of 
structure and physical properties (electrical conductivity, heat capacity, etc.) 
occur H]. Actually, the pressure - temperature phase diagrams of these sub- 
stances demonstrate two kinds of phase transitions. The first one corresponds 
to the sharp change of local order in the system. In solid state, these transitions 
are known as polymorphous ones. The possibility of polymorphous transitions 
in liquid area is under permanent discussion. Another kind of transitions, i.e. 
melting, is associated with global order changes. To describe the phase dia- 
gram, one has to describe both transitions in terms of one and the same sta- 
tistical model or, in other words, one has to describe all phases of the system 
in terms of one and the same set of order parameters. Such a description may 
be realized via local state representation in statistical mechanics of condensed 
matter introduced by Patashinski |^ . The idea that local state representation 
may describe the liquid - liquid phase transition was suggested in and we 
tried to realize it in the present paper. Here, we give brief introduction into 
Patashinski model (section 2), then we perform calculations for the system with 
two competing local structures. The results are applied to the binary system 
and to the system with two polymorphous local modifications, such as Se and 
S (section 3). Further applications of the model are discussed in Conclusion. 
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2 Patashinski model. 

Physical background of the model arises from the fact that in condensed sub- 
stance the interaction on atomic length-scale is strong enough to result in a 
substantial restriction of relative positions of particles. This is being referred to 
as the local order in condensed systems. The size of strongly correlated cluster 
and the nature of local order depends on the temperature, pressure, and com- 
position of the material. The local ordering implies that atoms positions in the 
local cluster correspond to some deformed ideal pattern. In close-packed mate- 
rials, one supposes several geometrical objects (the fragments of HCP and FCC 
lattices and icosahedron) to be such a pattern. The pattern can be character- 
ized by geometric parameters which are invariant under rotations (coordination 
lengths and numbers, characteristic angles in the pattern). These parameters 
are the same for most clusters and fix up the type of local order in the material. 

To account the limitations imposed on relative positions of particles by the 
local order, one describes the system in terms of local order parameters. Since 
the concept of local order assumes a space scale larger than mean interatomic 
distance, a coarsening of the description is necessary. The resulting model is then 
a simplification that allows one to understand global ordering starting from the 
assumed local order and from the interaction of local orders in different clusters. 
The formalism considered was first offered by A. Patashinski and M.Chertkov in 

and below in this section we will follow this paper. 

As was discussed above, material is treated as a system of small clusters 
of equal size, each cluster being in one of a set of N local structural states. 
The same set of local states is assumed for all clusters; a number i = 1,N is 
assigned to each of the states of this set. A state of the whole material is then 
described by the field cri(r), where cri{r) is an iV - component vector defined 
by the rule: tTi(r) 1, if the state of the cluster placed at point r has number 
i, and <Ji{r) = in other case. To proceed with a statistical theory of the 
local and global order, one has to find the non-equilibrium free energy (effective 
Hamiltonian, |^) as a function of the configuration {(Tj(r)}. One expects the 
effective Hamiltonian H{a} to include independent contributions from different 
clusters, and interaction energy that depends on the states of more than one 
cluster: 



Here and below, the Einstein's summation rule is supposed. In (|l|), ai is internal 
energy of a cluster having state number i, and interaction energy is written 
as two-cluster interaction of general form. Free energy (|l|) defines the order 
parameter field 




(1) 




{^} 
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Here, the partition function Z{ai{r),T) is a sum over all local states. As can 
be seen from (||), the Wi(r) is the probability of state numbered by i to occur at 
point r, and 

N 

J2w,{r)^l. (3) 

i=l 

Partition function Z may be treated as generating functional for the order pa- 
rameter field w{r): 

To calculate Z , one uses Hubbard - Stratonovich transformation. Namely, let 
us suppose that the quadratic form in Hamiltonian is positively defined, and its 
kernel may be written in the form: 

M,j{r)^E,,J{r). (5) 

(As it will be shown below, this is not relevant for the mean field approximation) . 
Introducing real iV-component "conjugated" field '(/'z(^); one rewrites Z as 

X n»;.rrexp[i(i?,,V.(r)a,(r) +a,(r)a,(r))]. (6) 

Here, the integration D'ip goes over all components ipi of the field in each 
point r. In (^), the summation over local states can be easily done: 

^ = /n.;.^^.Mexp{-^}, 
F = \ Er.r' J-Hr - r')£;.,^.(r)^,(r') - 
-TErHE^eMUEMr)+c^^ir))}]■ (7) 

In this relation, ^pi{r) is a continuous variable, so one is able to apply well- 
known method of analysis - to find out the most probable configuration which 
minimizes the thermodynamic potential, then to investigate fiuctuations around 
it, etc. Here, we shall restrict ourselves with the first step, which is the mean 
field approximation (MFA). 

Function ip{r), which minimizes the thermodynamic potential F in (0) (most 
probable configuration) obeys following equation: 

Er' J^'i^ - r')Mr') = Z-' exp[^iE,j^j{r) + «,)], 

Zo - E^ exp[^(S,,7/.,(r) + a,)]. (8) 

Substituting solution of this equation into (^, and using relation F{ip) — 
—T In Z, one gets the generating functional Z. The probability Wi (r) — —SF/Sai (r) 
coincides then with the right side of (||). Then, for Wi{r), one has 

w^ir) a^{r) J-i(r-r') < > ■ (9) 



3 



The F is similar to nonequilibrium thermodynamic potential in the Landau 
theory, with order parameter ip°'{r). At equilibrium, F is minimal. Let us 
define Wi(r) for nonequilibrium ip using (|^). Now the nonequilibrium potential 
F may be written in terms of Wi{r): 

-7^E.lnE,exp(i(E.' M,,(r - r')w,{T') + a,))]. (10) 

As can be seen from (^0|), condition (|^) is not essential in MFA. The nonequilib- 
rium thermodynamic potential written in terms of order parameter field u'i(r), 
is the base of Landau-like theory. At equilibrium, one has 

X (11) 



with normalizing condition (^. In explicit form, this equation may be written 



as 



u,,{r)^Z-\r)eM^) 

r' 

Z{r) = Y.eM^)- (12) 

1=1 

Here, the Ei{r) may be understood as the energy of cluster having local state 
number i, while the states of surrounding clusters are characterized by the mean 
probability Wj{r'). If the field a depends no on spatial coordinates, a{r) — Q!(G), 
then Ei,Wi are spatially homogeneous too: 

E, = M,j wj -f a, , M,j = M„ (r) . (13) 

r 

The normalizing condition (^ is fulfilled for solutions of (^2|) automatically. 
Depending on parameters of the theory (My , cti , T) , equations may have 
a set of solutions. In the vicinity of the first order phase transition, one ex- 
pects at least two minima of F, corresponding to stable and metastable phases. 
The lowest minima corresponds to the stable one. As an example of applica- 
tion of formalism described, let us consider N - local states model with special 
interaction: 

Mij = J5ij , 

a,(r)=0 (14) 

which is the Potts model In 0|, the Potts model was applied to describe 
the melting of a single-component material. The set of local states in this case 
should be understood as a set of allowed orientations of local clusters. The 
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physical picture is that each cluster of condensed material may be treated as a 
deformed cluster of ideal crystal. It is obvious for crystal phase; For liquid, this 
statement is an assumption which seems to be correct for a wide range of 
materials. For locally crystalline material, one studies the statistics of orienta- 
tions of local anisotropy. To obtain the effective Hamiltonian of orientational 
degrees of freedom, one uses following considerations. If two clusters share a 
boundary, then the minimum of their interaction energy corresponds to the co- 
inciding of their orientations. If the minimum is deep enough, then one is able 
to divide the space of orientations in equal cells and to use a coarsened descrip- 
tion of orientations in terms of orientation cells. The number of cells has to be 
chosen as to comprise the orientation attraction region, where the interaction 
energy is close to its minimum. Denoting the depth of the minimum as J, and 
the number of cells as N, one arrives at the effective Hamiltonian (|lj) of the 
Potts model which behavior is known. At low temperatures, the material is 
orientationally ordered: all clusters with probability « 1 have one and the 
same orientation (crystal). At high temperatures, all orientations are of equal 
probability (liquid). The melting of the material is therefore described as an 
orientation disordering. More detailed microscopic consideration of the relation 
between the Potts model and melting may be found in [|| . Note that an appro- 
priate choice of two phenomenological parameters J, N makes it possible do 
relate correctly several quantities, such as melting temperature, spinodal tem- 
peratures, the hidden heat of transition, heat capacity jump. That is more than 
one would expect from such a simple model. 

The melting theory, based on the Potts model, and our further considerations 
essentially deal with the statement that the local order does not change via 
melting - the same set of ideal patterns may be assigned to the liquid and to 
the crystal. This statement has no universal character, this is suggestion. In 
some cases, inherent amorphous structures may be involved during melting Q. 
But for systems considered, there exist supportive arguments for the crystalline 
local structure both in the crystal and in the liquid [||. 



3 System with two competing local structures 

In this section, we discuss an idealized model of a material having two possible 
local structures. Each structure has its own set of allowed orientations. One 
introduces two indexes to enumerate the local states: 

al; i = l,2; i = l^fc = l,7i; i = 2^k = l,m. (15) 

The upper index corresponds to the type of local order, while the lower one 
enumerates the orientations. The first and the second structures have n and m 
allowed orientations respectively. If only two-cluster interactions are accounted, 
then Hamiltonian takes the form 

n 

- ^ = " E E (^) + E ^fe (^)^^ki - r'K (/) , (16) 

r 2—1 r,r' 
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where a is the difference of internal energies of two local structures. The kernel 
of interaction differs from zero only for nearest neighbors. We suppose it may 
be written in the form: 



Mll = Ji6ki; Mit = J25kV, Mlf^Mf,'=e (17) 

Such a form of interaction arises from following. Let the two structures dif- 
fer sharply from each other. If two clusters with different structures share a 
boundary, then their interaction energy is e despite of their orientations; the 
orientational interaction of neighbors having one and the same local structure 
is described by the Potts model. Parameters Ji, J2 are the depths of orienta- 
tional interaction of structures 1 and 2 respectively, and e - is an interstructural 
surface energy. Following notations will be used below: 

e = ev; Ji = Jiv; J2 = J2J^, (18) 

were v is the number of nearest neighbors. 

For the mean - field order parameter, one supposes, in analogue with the 
Potts model: 

1 1 p-wi 

^ n — 1 

2 2 , , 

wi W2; Wfe^i = — • (19) 

Here, p is the mean probability of the first structure. Relations ( [l9|) realize 
the idea of orientational ordering; the number 1 is assigned to the orientation 
which is most probable in the crystalline state. Now one is able to introduce 
following classification of possible phases in the model: 

1. Wi ^ p, p 1 - crystal with the first type of local structure; 

2. wi = p/n, p ^ 1 - liquid with the first type of local structure; 

Z. W2 ^ ^ — Pt p ^ Q - crystal with the second type of local structure; 
4. u'2 = (1 — p)/m, p ^ - liquid with second type of local structure; 
For the thermodynamic potential ( p^ ) per one cluster, one gets 



f= - J2{1~pY\W2 



2 2(m-l) 
(l-^i)' 



r^-T 2(n-l) 
ep{\-p)+T{\~p) ln(l -p) +rplnp 
Tplnwi + T[l — p)\ti'W2 + pa , (20) 
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where the vahies wi —wi/p, W2 — ^2/(1 — p) obey equations 



Wi 



W2 



1 + (n — 1) exp 
1 + (m — 1) exp 



Jip(l — nwi) 
T{n - 1) 

T{m-1) 



(21) 
(22) 



With the help of (po|), the phase diagram of the model can be plotted in MFA 
for two physically different cases, considered below. 

The first case is binary system in which the second type of local order occurs 
due to presence of second component, i.e the binary system with limited misci- 
bility of the components. The two types of local order correspond to the pure 
substances A and B. The model can be applied to systems in which the two 
types of local order differ sharply (have different groups of local symmetry or 
incompatible interatomic distances). In detail, the calculations are performed 
m 0, where it was shown that all types of phase diagrams for binary sys- 
tems with limited miscibility of the components (eutectic and monotcctic ones) 
are described correctly. To demonstrate validity of the model, we calculated 
the concentration - temperature phase diagrams for well known binary systems 
Ga-Pb and Ag-Cu. The results are plotted on fig.l. One can see some devia- 
tions from experimental data published by Elliott [Tl|] . Determination of model 
parameters was made by following way. The parameters are Ji,n, J2,m,a,e. 
Parameters Ji , n and J2 , m correspond to pure components and were deter- 
mined from the heat capacity temperature dependencies of pure Ag, Cu, Ga, 
Pb [ p^ . For binary system, where the concentration of second component is an 
external parameter, the value of a does not play any role: it should be chosen as 
to provide the fixed concentration of mixture. The only parameter to fit known 
diagram is e! Thus, the agreement with the experiment is rather good. 

For the Ag-Cu system the parameters are : Ji = 1.356, n = 4 (Cu ), 
J2 = 1.234, n = 4 (Ag), e — —0.275. The Ag melting temperature was 
chosen as the unity. For Ga-Pb, the parameters are Ji — 3.18, m — 11 (Pb), 
J2 = 1.6, n = 11 (Ga), e = —1.596. Temperature unity corresponds to the Ga 
melting point. 



4 Selenium and sulfur. 

Consider now the system where atoms rearrangement changing the type of local 
order in small cluster is allowed. Instead of binary system, one has the substance 
with polymorphous modifications. The thermodynamic potential (^0|) may be 
applied to plot the phase diagram also. Phase transitions arise, if one changes 
the model parameters - Ji, J2, a, e. In physical situation, the change of param- 
eters may occur under pressure. Here, we describe the pressure - temperature 
phase diagrams of Se and S which are known experimentally |^ . The coefficients 
in the model were chosen as to fit the experiment. The calculations are presented 
at fig. 2, fig. 3. For sulfur, the calculations were done by following way. First, the 



7 



parameter n (number of orientations of low - pressure structure) was determined 
from the heat capacity data for melting at atmospheric pressure. Since there is 
no reliable heat capacity data for high pressure structure, the parameter m was 
supposed to be near 4. We suggest this value because of our experience in simi- 
lar calculations: for systems with high density, the number of orientations tends 
to two "magic" values, 4 or 11. In our case, the value 4 provides better corre- 
spondence to experiment. We suppose linear pressure dependence for Ji,J2,a 
and no pressure dependence for other parameters ( n,m,e ). The coefficients in 
linear dependencies Ji = aip + bi, J2 = a,2P + 62, a = 03^ + 63 were determined 
from the "tales" of experimental diagram: ai, bi - from the low temperature 
and low pressure line between phases 1 and 2, 03, 63 - from the interphase line 
near critical point, 02, 62 - from the melting line at high pressure. The central 
part of the diagram was then plotted automatically. For sulfur, the parameters 
are n = 2.5, to = 4, ai = 0.527, bi = 2.433, 02 = 0.558, 62 = 3.296, 03 = 
—0.094, bs — 1.535 ,e = —4.506. One can see good correspondence between 
calculations and experimental data. At higher temperatures, there exist another 
phase transition into metallic liquid, but it is out of our consideration. 

For selenium, the calculations were done by similar way. The only difference 
is that the position of melting line at high pressure is not known from exper- 
iment, so we tried to reproduce the interphase line between crystal phases at 
low temperatures, which position on the diagram was supposed by V.Brazhkin 
So, the correspondence to the experiment in the case of Se allows to wish 
better. The parameters for selenium are n = 2.5, m = 4, ai — 0.623, bi = 
2.433, a2 = 0.443, 62 = 1.3, a3 = -0.226, 63 = 1.145 ,e = -4.133. 

In both cases, temperature scale was chosen as to make the melting temper- 
ature at atmospheric pressure equal to unity. 

5 Conclusion 

The main advantage of the model considered is the fact that all possible phases 
are described by one and the same order parameter - vector of mean probabil- 
ities. This gives one possibility to write out the expression for thermodynamic 
potential which is valid for all situations, and thus to investigate the phase 
transitions. 

In the model, the drastic changes of structure in liquid are described as the 
continuation of the line of polymorphous phase transition into the liquid area, 
where it terminates in the critical point. Below critical temperature, the two 
types of local order are well distinguished, and phases 2 and 4 (see classifica- 
tion) are divided by the line of the first order phase transitions. At critical 
temperature, the transition is of second order. Above it, the local structures are 
degenerated by thermal fluctuations. Note, that the well-known phase diagrams 
at fig.l and phase diagrams presented at fig. 2, fig. 3 are the diagrams of one and 
the same model. The only difference is that for fig. 2, fig. 3 one uses pressure 
instead of concentration as a thermodynamic parameter. 

Further application of the model concern an interesting line corresponding to 
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the equilibrium of liquid and crystal having different types of local order (phases 
1 and 4) . We hope the model to explain " anomalous" melting of covalent crystals 
with diamond structure (C,Si,Ge) and compositions of type AjuBy (In-Sb, 
Ga-Sb) with zincblende structure. Melting of these substances is accompanied 
with changing of the type of chemical bonds and sharp changing of local order 
(coordination number rises approximately in two times) jl3| This results 
in essential increasing of density and in the semiconductor - metal transition 
during melting At the same time, at high pressures, these substances 

transform into metallic crystals with the structure of white Sn One can 
say that there exist a line of semiconductor - metal transition which changes its 
direction at the crossing point with the line of metal crystal melting. This point 
is the point of three - phase equilibrium of covalent crystal, metallic crystal, and 
metallic liquid. In principle, the line of coexistence of crystal and liquid with 
different local structures may terminate at lower pressure in the crossing point 
with the line of covalent crystal - covalent liquid transition. Then the line of 
polymorphous transition should be plotted into the liquid area and terminated 
in the critical point. For substances mentioned, this is doubtful because the 
"covalent" melting should take place at large temperatures It is more 

reliable that the line terminates near the liquid - gas equilibrium, which is 
out of the validity of the model due to the destroying of local order at high 
temperatures. 
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Figure 2: Pressure - temperature phase diagram of sulfur. Solid circles represent 
the experimental data, solid lines correspond to our calculation. Critical point 
is shown as an empty circle. The enumeration of phases corresponds to the 
classification given in the text. 
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Figure 3: Pressure - temperature phase diagram of selenium. Solid circles rep- 
resent the experimental data, solid lines correspond to our calculation. Critical 
point is shown as an empty circle. The enumeration of phases corresponds to 
the classification given in the text. 
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